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The small perturbations method has been extensively used for waves scattering by 
rough surfaces. The standard method developped by Rice is difhcult to apply when we 
consider second and third order of scattered fields as a function of the surface height. 
Calculations can be greatly simplified with the use of reduced Rayleigh equations, 
because one of the unknown fields can be eliminated. We derive a new set of four 
C ' reduced equations for the scattering amplitudes, which are applied to the cases of 

a rough conducting surface, and to a slab where one of the boundary is a rough 
surface. As in the one-dimensional case, numerical simulations show the appearance 
of enhanced backscattering for these structures. 
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I. INTRODUCTION 
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O^ , The scattering of electromagnetic waves from a rough surface has been studied in 

Cn such as radio-physics, geophysical remote sensing, ocean acoustics, and surface optics.H 
^Q One of the earlier theory used, is the-small perturbations method (SPM) originally developped by Rice,ll^ 
^C. this theory still remains of interestllilllj'tJ because perturbative terms of order higher than one can pro- 
^— V '■ duce enhanced backscattering, or improve predictions accuracy in an emission model. Although, the Rice 
"T*-, method can be used in principle to determine all orders in the perturbative development, very few works 
C^ use terms of order higher than one for a two-dimensional surface due to the calculations complexity. The 
C ', second order has been written in a compact form by Voronojrichllj in his work on small-slope approxi- 
I ' mation, and only recently the third order has been presentedJ13 However, there exists a different way to 
' O obtain (SPM) which dates back from the work of Brown et alt3 Using both the Rayleigh hypothesis and 
5 , the extinction theorem, they have obtained an integral equation, called the reduced Rayleigh equation, 
which only involves the incident and the scattered field alone. In their method the field transmitted 
through the surface has been eliminated in such a way that the scattered field becomes a function of 
_ j^ ' the incident field only. This reduced equation.. has, been extensively used by Maradudiri-ei al. to study 
k> , localization effects by a conducting surface ,EIlE3'E3 coherent-effects in reflection factor,E] scattering by 
5_] ■ one-dimensionalEil and two-dimensional conducting surface.E3 It has to be noticed that the third order 
Ci ■ perturbation term was already explicited in the work of Ref. G2. . .. . 

In recent years, similar studies have been done in the case of thin films bounded by a rough surf ace, llirE3 
but only in the one-dimensional case.^^ In order to calculate the two-dimensional case it becomes neces- 
sary to derive an extension of the reduced Rayleigh equation for this system. In the present paper, we 
first study a surface where a down-going and up-going fields exist both on the upper and the bottom 
side of the film. We show the existence of four equations, that we also call reduced Rayleigh equation, 
they have the property that one of the down-going or up-going-^elds has been eliminated. With these 
equations, we rediscover the equation obtained by Brown et alrB for one rough surface, and derive the 
corresponding ones for a slab where one of the boundary is a rough surface. Next, a perturbative devel- 
opment up to third order is obtained in a compact matrix form for these two systems. This third order 
term is mandatory if we need the expression of the cross section up to fourth order approximation. As 
in the one-dimensional case, the results of the incoherent cross section show a well defined peak in the 
retroreflection direction. 

In the case of s rnall p r pughness metallic surfaces ttis peak was originaly explained by the infinite per- 
turbation theory,tZHl§ll3 and further developments^ have shown that the major contribution for the 
enhanced backscattering peak comes from the second order term in the field perturbation. However, 
in the one dimensional case the enhanced backscattering for a rough surface, appears only for a [TM) 
incident wave due to the fact that plasmons polaritons only exist for this polarization. In the two dimen- 
sional case, due to the existence of cross-polarization, we will show that an incident {TE) wave can excite 
a {TM) plasmon mode which can transform into a {TE) or {TM) volume electromagnetic wave. Thus, 



the enhanced backscattering is present independently of the polarization of the incident and scattered 
waves. For a rough dielectric film bounded by a conducting plane the enhanced backscattering is present 
for both (TE) and [TM] incident waves, even in the one-dimensional case because guided waves exist 
for these two polarizations. The qualitative effect of the two-dimensional surface,is,particularly sensitive 
when we study thin film, for instance, in one-dimensional case, satellite peaksHilS appear on each side 
of the enhanced backscattering peaks, however, in the two-dimensional case, the coupling between {TE) 
and [TM) modes attenuates drastically these peaks. 

The paper is organized as follows. In Sec. EB, we derive the four new reduced Rayleigh equations. In 



Sec. flj, we introduce the diffusion matrix, and in Sec. IV, we determine the perturbative development 
up to third order term in the surface height, for a rough surface alone, and for slab with a rough 
surface located at one of the boundaries. In Sec. ^, we introduce the Mueller matrix and the definition 
of the statistical parameters for the rough surface. We then obtain the bistatic matrix in terms of a 
perturbative development. Numerical examples which show the enhanced backscatteri ng a re presented 



in Sec. VI. Conclusions drawn from the results of our calculations are discussed in Sec. VII 



II. DERIVATION OF THE REDUCED RAYLEIGH EQUATION 

The reduced Rayleigh equation was obtained for a two-dimensional surface by Brown et alX^ using the 
extinction theorem and the Rayleigh hypothesis, it allows to calculate the scattered field from the rough 
surface. Now, if we want to compute the transmitted fpid from the rough surface, one has to introduce 
an other reduced Rayleigh equation derived by Greffet.EJ However, these two equations were established 
in the case where there is no up-going field inside the medium, thus they cannot be used to obtain the 
field scattered by a slab with a rough surface in its upper side. In fact, to generalize these equations to 
a slab, we have to consider all the fields shown in Fig. ^ We will prove that there exists four reduced 
Rayleigh equations, which involve only three of the participating fields Eq, Eq, Ei , Ef. 

We consider that each electromagnetic waves propagates with a frequency w, and in the following the 
factor exp{~iujt) will be omitted. We choose to work with a Cartesian coordinate system r — {x, z) — 
{x, y, z), where the z axis directed upward, and we consider a boundary of the form z = h{x). Moreover, 
we suppose that there exists a length L for which h{x,y) = 0, if \x\ > L/2 or |y| > L/2, L may be 
arbitrary large but finite. 

A Propagation equations and boundary conditions 

The electric field E satisfies the Helmohltz equation in the two media: 

(V^ + eoK^)E"{r) = for z > h{x) , (1) 

(V^ + eiK^)E\r) ^ ior z < h{x) , (2) 

where Kq — lo/c. Since the system is homogeneous in the x — {x,y) directions, we can represent the 
electric field by its Fourier trapsform. Thus, using the Helmohltz equation, we deduce the following 
expression for the electric fieldclllj in the medium : 

E"ir) = I ^E'-{p) exp(ife°- ■ r) + J ^E'+ip) exp(ifc°+ • r) , (3) 

where (see Fig. ^) 

aoip) ={eoK^ - p')K (4) 

kp"^ =p ± ao{p)ez . (5) 

In fact, when writing such a definition, we made implicitely the assumption that the Rayleigh hypothesis 
is correct. This representation is only valid when z > max[h{x)] and in that case E'^~ (p) represents the 
incident wave amplitude. In order to be correct we need to add an explicit dependence in the z coordinate 
like (see Ref. ||): 

E°- =E°~{p,z) £;°+ = £;°+ (p , z) (6) 



However, explicit calculations in the case of infinite conducting surfaces]^ and for a dielectric mediuma 
without this hypothesis, have shown that the perturbative developments are identical. Tht^alidity of 
this hypothesis is by no doubt a matter of convergence domain as discussed by Voronovich.oEj 
In the medium 1, we have a similar expression : 



^'^"■^ = / ^^'^^^^ "^P^^^^" ■ ""^ + / (^^'^^^) """P^^^^^ ■ ""^ ' (^^ 



fei±=p±ai(p)e,. (9) 



^^^^-(p)exp(ifc-.r) + /|^^, 

where 

Qi(p)=(eiX2-p')^ (8) 



We decompose the vectors E{p) on a two-dimensional basis due to the fact that V.E{r) = 0, which gives 
the conditions : 

fc°±.£;°±(p) = 0, kl^.E'^{p) = Q. (10) 

Then, we define the horizontal polarization vectors H for (TE) and V for (TM) in medium by : 



eg(p)xfc°^ _^ ao(p) . IIpI 



\\eH(Po 
with similar expressions for the medium 1: 



'""''" life xfc^ll '"'"^- "" 

4V).#^^^^^^±^P-^^.. (14) 



|eH(Po)xfepll x/el^o" V^i^o 



p 



So, we decompose the waves in medium on the basis [p]° = (ey (p),eH(p)), and [p]'^^ = 
(4+(p),6hW): 



V ^H (P) / [p]0- \ ^H \P) J [p]0+ 

and for medium 1 on the basis [pY~ = [e-^ [p), enip)) and [p]^^ = (e^ (p), enip))'- 

The electric E{x, z) and magnetic fields B[x, z) = 4-V x E{x^ z), satisfy the following boundary condi- 
tions : 

n{x) X [E°{x, h{x)) - E^{x, h{x))] = , (17) 

n{x) ■ [eoE°{x, h{x)) - eiE\x, h{x))] = , (18) 

n{x) X [B°{x, h{x)) - B\x, h{x))] = , (19) 
n{x) = e^ — Wh{x) . 

Let us introduce the fields Fourier transform, Eq. (gj) and Eq. (0), in the boundary conditions Eqs. ( |r7H19| ), 



they give : 



E / (^"(^) >< E'^iP)eM^K^■r.) = E / (^"(^) >< E"^{p)eM^kl; ■ r^) , (20) 



a=± 



a=±' 



£0 

e 



; E / (^"(^) ■ ^"''(p) ^^p^^'^p" • ^-) = E / (^"(^) ■ ^'"(p) ^^p(*^p" • ^-) ' (21) 



E / (|^"(^) ^ K X ^°"(p)] exp(*fcp" • ^-) = E / (|^"(^) X t^p" X ^'"(p)] ^^p(*^p" • ^-) ' 



r^ = a; + /i(a;)e^ , kp"" = p + aao{p)e^, k^ = p + aai(p)e^ , 



(22) 
(23) 



where the summation includes the two possible signs: a = ±, hnkcd to the propagation directions. We 
will also use the condition V • E'^{x, z) = V • E^{x, z), which gives the relation 

E / ^K'' ■ ^""(P) exp(*fc^° ■ ^-) - E / ^''p ■ ^'"(P) ^MK"" ■ ^-) • (24) 



B Fields elimination 

The equations (p0|-p2|) and ( p4|) , are all linear in the fields £'°^, E^^ , E^^ , E^^ . In order to eliminate 
E^~ or E^~^ in the equations ( pO[p^ ) and (|2^), we will take the following linear combination of their left 
and right members: 



Jcfx [kl' X (Eq. {^) + {Eq. (|2|)) - k^^{Eq. (^) - n{x)(Eq. (H))] cxp(-zfei^ • r. 



(25) 



with fc^ = M + bai{u)ez, and where b = ±, has to be fixed according to the choice of the field we want 
to eliminate. With the vectorial identity, a yi {b 'k c) ^ b{a ■ c) — c{a ■ b), the right member of Eq. ( psf ) 
can be written : 



E 

a=± 



.2 d^P 
d X- 



(2^)2 



K' + kl^) ■ nix) E'^p) + (kl' kl^) ■ E'^p) nix) 
-nix) ■ E'^ip) ikl' ~ kl^)] eM-^(.kl' - kl^) ■ r.) . 



(26) 



We have now to discuss the different cases depending on the relative sign between a and h: 

1) If Q = —6 , we can use an integration by parts (see Appendix ^) to evaluate n(a;) = e^ — Vft.(x). 
Then we can make the replacement : 



nix) ~ e-z — Vhix) < — > n(a;) = e^ 



iu-p) 



ihaiiu) — aaiip)) 



(27) 



It has to be noticed that the denominator ihaiju) — aaiip)) does not present any singularity because 
a = —b. For the first term in the integral (|26| ) we obtain: 



(fc- + fe-).n(a.)£;-(p) 



-bE^^ip) 



iaiiu) + aiip)) 
0. 



[u^ -p^ + aiiu)^ - aiipY' 



(28) 



The last equality can be easily checked using Eq. (||). For the sum of the second and third terms of 
Eq. (p6|), we have also: 



(C - fcp") • E^'iP) nix) - nix) ■ E'%p) (fe^^ - k^^) = , 



(29) 



due to the fact that: 



n{x) = 






bai{u) — aai{p) 



(30) 



2) If a — b , we can use agam the integration by parts only if ai{u) ^ ai{p). Then we have to consider 
three cases: 

2-a) u j^ p and u ^ p, as in the previous case by using an integration by parts we show that Eq. ( P6[ ) 
is zero. 

2-b) u=p, thenk]^ = k}^: 

- f d'xikl' + fc^") • n{x) £;i"(p) exp(-i(fc^'' - fc^") ■ r^) ^ ~ j A^x2k'^ ■ n{x) E^%p) 

= —2hai{u) I dx 
= -2hai{u)L^ , 
because J d'^x\7h{x) — 0, and 

{kl' - kl^) ■ E'%p) n{x) - n{x) ■ E'^{p) (fc^" - k^;) ^ . 
2-c) u = -p^Q, then k]j; - fc^" = 2m , 

/ d^x n{x) exp(-i(fc" - fc^") • r^;) = / d^x (e^ - Wh{x)) exp(-2m • x) 

= Bz d^xcxp{—2iu ■ x)x — Bx / dy [exp(— 2m • a;)ft,(a;, y)]^~_'^ ,2 

-By dx [exp(-2m • x)h{x,y)]lZ^^]^/2 

— e.z5{u) when L -^ +oo 

= since it 7^ . (33) 



(31) 
(32) 



This result implies that the expression (|26|) is also zero in that case. 
We can summarize all the above results in the form : 

^Jd'xikl^ + kl^) ■ n{x)E^%p) exp{-i{kl' - k^^) ■ r^) = -2bai{u) Sa.bS^^pL^ 



= -2bai{u) Sa,b (27r)2 d{u - p)E^''{u) (34) 
when L — > +cx) , 



where 6u^p is the kronecker symbol, and S{u — p) = (27r)^/i^ 6u,p the Dirac function 
After an integration on p and a summation on a, we obtain for the expression (26) 



lb/ 



-2bai{u)E'-''{u). 

We see that we can eliminate the field E (m) or E ^(u) depending on the choice made for & = ±. 
Now, if we consider the left member of Eq. (Ea), we have: 



(35) 



E 

a=± 



d^a; 



d^p 

(2^ 



{kl' + kl^) ■ n{x) E'^p) + {kl' kl^) ■ E'^p) n{x) 

exp{-z{kl' - fe;^'') . r^) . 



-n(a:).£;-(p)(-^fc--fc-) 



Using an integration by parts, we replace n{x) by (27) 

(u-p) 



n{x) 



UU ,0a 



{bai{u) - aaoip)) ibai{u) ~ aao{p)) 



(36) 



(37) 



In this case there is no need to discuss the relative sign between a and b because bai{u) — aao{p) ^ 0, 
due to the fact that eo ^ ci- We then obtain : 



Oaf 



{K + K)-n{x)E"%p) = - 



bai{u) — aao(p) 



where we have used the definitions fl) and (0) . The remaining terms of Eq. 



E'>%p) 



I give 



(38) 



{K 



kl^) ■ E^^ip) n{x) n{x) ■ E"^{p) ( Jfe^" - fc°'') 



(fc 






p 



£;"'^(p) n(a;) - n{x) ■ E^^p) (fe" - k^^) + n{x) ■ E^'%p) [k 



.16 



ei 



bai{u) — aQ;o(p) 
Introducing the following notation : 



(ei 



121 k^b 



£1 



/(a|p) = / d^x exp(— ij> ■ x — ia h{x)) . 



(39) 



(40) 



and taking into account the expressions (p5|), 
in the form: 



i9|), we express the resulting linear combination (^ 



E 



d^p /(6q!i(m) — aao(p)|M — p) 



(2^)2 



bai{u) ~ aaa{p) 



if^£;°''(p)--^(fe 



Ifc _ . OaN 
« "'p ^ 



£;"°(p) 



26ai(M) 
(ei - eo) 



E'^iu) 



(41) 



this expression represents in fact a set of two equations, depending on the choice for 6 = ±. The last step 
is to project (O) on the natural basis of E (it), namely [m]^'' = {ey (it), eniu)), which has the property 
to be orthogonal to fe^ , so it eliminates the second term of (B^) l.h.s. Let us notice that in order to 
decompose E "(p) on [p]"°, one has to define a matrix M ' {u\p) transforming a vector expressed on 
the basis [p]°" into a vector on the basis [u]^'', multiplied by a numerical factor (eoei)^ Kq introduced 
for a matter of convenience: 



\eH{u)-ey{p) eH{u)-eHip). 



(42) 



C The reduced Rayleigh equations 

With the definitions (O-nJ) the matrix M takes the form: 



M ' {u\p) 



||m||||p|| + abai{u) ao{p) ii ■ p —be^ Ko ai{u) {u x p) 
*^ ae^ Koaoip){ux p)^ {eoei)^ K^ii-p 

and, the two reduced Raleigh equations resulting from Eq. (EI|) read : 



a=± 



d^p I{bai{u)-aao{p)\u~p)—-ib,Oa 



(2^)2 



bai{u) — aao{p) 



M ' {u\p)E"''{p) 



(fi - eo) 



u), 



(43) 



(44) 



where we suppose that E^°'{p), E^ (u) are respectively decomposed on the basis [p]*^" and [u]^''. We can 
derive a similar equation where E is now eliminated, by simply exchanging cq and ei in ( [43|) and ([44|), 
due to the symetry of the equations (0J^),(O|l9), we get : 



(45) 



d^p I{baoiu) ^ aai{p)\u - p) ^^bsa, , . „ia. x 2 b (eg e^) 2 ao{u) n^ 

^^ t 'PTyi 1 — r~\ r~\ ^ {u\p)E (p) = E (m), 

^^ I {2TTy bao{u) - aai{p) (ei - eo) 



a=± 



'M°'''^" (u\p) = [ ll"IIIIPJI+"^"o(«)ai(p)w-P -bel Koao{u){uxp)A 
\ ael Koai{p){uy.p)z {eoei)^ K^u-p J 

In the next sections, we will show how these equations greatly simplify the perturbative calculation of 
plane waves scattering by a rough surface. In order to keep a compact notation, we introduce new 
matrices M^ given by : 

—-ibfia, , , I(bai(u) ~ aao(p)\u — p)-—ib.Oa, , , ,,„, 

Mh {up)= ' ' "^^;' ^^ M ■ MP , 47 

Af^'''\u\p) ^ I^ba,[u)~aa^{p)\u^p)^bM^^^^^ ^^g^ 

hao{u) - aai{p) 

III. THE DIFFUSION MATRIX 

We are interested by the diffusion of an incident plane wave by a rough surface from the previous 
formalism. We define an incident plane wave of wave vector A; ~ as: 

^ Pa 

E^-(p)^{2^f5{p-p^)E\p,). (49) 

We are naturally led to introduce the diffusion operator R: 

E'+{p)^R{p\p,)-E'{p,), (50) 

which can be represented in a matrix form, using the vectorial basis described above : 

Kfr.l'n \ - fRvv{p\Po) Rvh{p\Pq) 

R[P\P0) - [R^y(p\p^) Rhh{p\Po)j [,-j^[p+] 



The field in medium is now written (using the decomposition (p^ 

E\r) = E\p,) exp(ifc°; ■ r) + J ^R{p\p,) . E^{p,) exp(ifc°+ • r) . (51) 

IV. A PERTURBATIVE DEVELOPMENT 

In order to obtain a perturbative development, one has to make a perturbative analysis of the given 
boundary-problem. A direct approach which uses an exact integral equation named the extended bound- 
ary condition (EBC) (see Ref. [4|,p6|) requires tedious calculations. An other issue is to us q thC pB-ayleigh 
hypothesis in the boundary conditions. This is the method generally used to obtain (SPM)ll3ll3'Ej. But a 
great deal of simplifications_ean be achieved if we are only interested by the field outside the slab. It was 
discovered by Brown et alr3, that an exact integral equation can obtained (excepted for the Rayleigh 
hypothesis), which only involves the scattering matrix R{p\pq). The proof is based on the extinction 
theorem which decouples the fields inside and outside the media. In this section, we will show how to 
obtain this integral equation from the previous development, including a generalization to the case of 
bounded random media. We seek for a perturbative development of R in power of the height h: 

R{p\Po) ^R^"\p\Po) + R^'\p\Po) + R^'\p\Po) + R^'\p\Po) + ■ ■ ■ ■ (52) 

One can easily prove, that this development takes the following form (see Appendix pf): 

RiplPo) - {27rf6{p^Po)X^'\po) + aoiPo)x'''\p\p,)hip-p,) 

/r\ ly ("P^ 
7^^ iP\Pi\Po) Hp - PiMPi - Po) 

+ao(Po) // -n^J^^^ {p\Pi\P2\Po) Hp ~ Pi)HPi - P2)HP2 - Po) , (53) 



where h{p) is the Fourier transform^'^ of h{x): 

Mp)^/c1^.cxpHp..)M.). 



(54) 



We will now exemplify the power of the reduced Rayleigh equation for the three configurations mentioned 
in the introduction. 



A A rough surface separating two different media 

We consider a rough surface delimiting two media which are semi infinite, see Fig. 0. We suppose 
that there is no upward field propagating in the medium 1, so E ^ =0. With the choice, 6 = +, in 
equation (44), we obtain the following integral equation for the scattering matrix Rg ^^^^-^^pIpq) for a 
single surface'^^ (the subscript s means a single surface located at z = 0) 



d^P -ji-rl + '" 



i+,o- 



(27r) 



M^ ' {u\p)-Rseo,eAP\Po)+Mh ' {u\po) ^ 



(55) 



Hit 



(This equation has been already obtained making use of the extinction theoremJl3 It has to be noticed 
that since the r.h.s of ( |55| ) is null, one can simplify the second lign of the matrices, M ' ,M ' , by 
a factor {ei)^ , then they coincide with the M,N matrices derived by Celli et a£3.) 

In order to construct a perturbative development, the method is simply to expand in Taylor series the 
term exp{ia h{x)) inside I{a\p), (Eq. (Ec 



^'^ ia^ 



I{a\p) = (27r)2 5{p) - za /i^) (p) - ^ /i^^) (p) _ !^ /,(3) (p) 
h^"\p)= f d^xexp{^ip-x)h''{x), 
and, to collect the terms of the same order in h{x). Let us define the matrix 






(56) 

(57) 

(58) 



the classical specular refiection coefficients for (TM) and {TE) waves are given by the diagonal elements 
of the matrix 



^ttIO 



y (Po) = -Dio(Po) • -Dio(Po) 



^seo,e^iPo) 



ai(Po) + ao(Po) 



-V'\P0), 



and 



(59) 



Introducing ( |5q ) in (pq), we obtain for R^^g^^-^ a perturbative development of the form (|53|), where the 
coefficients are given by: 



(60) 

(61) 
(62) 



^/eo,ei("IPilPo) = ai(w) Q^(m|Po) + MPo)Q (w|Po) - 2P(m|Pi) • Q^ (PiIPq) , 

^.eo,eiHPi\p2\Po) ^ --^ [(Oj (m) + Oq (pg)) Q ("IPo) + 2 ^1 (m) ftQ (Po) Q ("IPo) 

+ iPiu\p,)xfl^^Ap,\p2\po) + I {ai{u) - aoip^j) ■ PHp^) ■ ^''(PalPo) , (63) 



with 

Q (w|Po) = 5 — 7-1 [M {u\u)] ■ [M (m|po) ± M (m|po) • X (po) , (64) 

2ao(Po) 

P(m|Pi) ^ (ai(M) - aoiu)) [m'+'°+(«|«)]-i • m'+'°+(m|pi)] (65) 

and, after some simple algebra we obtain: 

Q+(«|Po)-(^i-eo)[^^oN]"'-f''"''""^^°"""'"'^"^"'^^"^'^'^° -6o'^o«i(«)(Hxp„)A 

■ [DtoiPoT' , (66) 

"o(Po) 

(<^oaiiPo)\\u\\\\po\\-eiai{u)al{po)u-po -^e§ Koai{u) ai{p„) {ii x Pq)A r n+ . ni^i 

V -eo ei^oao(Po)("xPo)^ ^o"i(Po)"-Po / 

(67) 

P(Tx|pi) = (ei-6o)[:D+(«)]-i. fll"lll|Pill + "i(")"°(Pi)^-^i -4^oai(«)(nxpi)A ^^^^ 

It can be easily checked that X is the well known first order term in perturbation theory which was 
obtained by Rice.ll3 AftejE_some lenghty calculations, we have proven that Eqs. ( p2| , |63[ ), are identical to 
those found by Johnson.Ej Thus our expressions (62), (|63|), are a compact manner to write the second 



and third order terms of the perturbative expansion, moreover, they are well adapted for numerical 
computations. However, it has to be noticed that only the first term X is reciprocal. Since the second 

(2) (3) 

and third-order perturbative terms are included in an integral, the coefficient X , X , are not unique, 
however they can be put into a reciprocal form (see Appendix |B| ). 

It is worth to notice that we can follow an analogous procedure to calculate the transmitted field. By 
taking 6 = — in Eq. (p5|), we get: 

/ p^M^r'-Hp) • E^-ip) = iifplii^ ^o-(.) . (69) 

J (27r)^ (ei - Co) 

This equation was already obtained with the extinction theorem.ci 

B A slab with a rough surface on the bottom side 

We consider a slab delimited on the upper side by a planar surface and on the bottom side by a rough 
surface, see Fig. y. Since there is no incident upward field in medium 2, the scattering matrix obtained 
in the preceeding section is sufficient to determine the scattering matrix of the present configuration. In 
order to get a proof, let us introduce some definitions as explained in Fig. ^. The scattering matrix for 
an incident plane wave coming from the medium 0, and scattered in the medium 1 is given by: 

v"{p\po) = i2nfSip-p,)v'\p,), (70) 



^10 



where V is defined by (p9|). The transmitted wave in the medium 1 is given by: 

T°(p|Po) = i2nr 5{p - p,) ^ t'\p,) , (71) 

ai(Po) 

t'\p,) ^ 2a,{p,) f(^«^i)^ 0) .[D+ (po)]-i . (72) 
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Now, when the incident wave is coming from the medium 1, we have similarly : 



10, 



V {p\po) = -(27:)' S{p-p,)V (po), 



(73) 



T ip\p,) = i27:fS{p-po)T'ip,). 



^H 



The scattering matrix R^ ^ ^ for the rough surface h which is located at z 
media of permittivity ei and €2 is given by: 



(74) 
H ^^ , and separating two 



R 



H 



(pIPo) =exp(i(ai(p)+ai(po))^)-Rs£i,e2(p|Po)^ 



(75) 



where the phase term comes from the translation z = —H, (see (B2)), and Rs ei,e2 denotes the scattering 
matrix Rseo,ei of the previous section, where we have replaced eq by ei, and ei by 62. Furthermore, if we 
define the product of two operator A and B by: 

d'Pi -IT, 



iA-B)ip\p,) 



{2ny 



■ A{p\Pi) ■ B{Pi\Po) . 



(76) 



one can easily prove for the configuration shown in Fig. Q that (we use for the fields the notations of 

Fig.g), 



E 



1+ 



<,.-^-^" 



Rf.,,2-V'-E'^. 



= V -E" 



T -E 



1+ 



(77) 
(78) 



where E^^ (p) ~ (Stt)^ S{p —Pq)E^{pq). Thesa-equations have been recently used to calculate in the first 
order the field scattered by a layered mediumJlil In fact, as we shall see below, these equations allow us 
to obtain all orders of the field perturbation. The expression (IT^) is analogous to the Dyson equation 
usually used in random media.u So, we are naturally led to introduce a scattering operator U : 



^n 



:flf,,,,,.c7.T".£;"- 



(79) 



which satisfies the equation 



U 



^'•«f.„e.-C/- 



(80) 



If we define by RdiplPo) the global scattering matrix for the upper planar surface and the bottom rough 
surface by. 



E^'^ip)=R,ip\p,)-E\p,)., 



with Eqs. (^8 79), the scattering matrix becomes: 



Rd = V 



T^<,.-I7-T°. 



(81) 



(82) 



We can improve the development (80), by summing all the specular reflexions inside the slab, this can 



<o) 



be done by introducing the operator U which satisfies the equation 



U 



(0) 



i+r-R^l-u'"' 



(83) 



where R^ ^ ^ is the zeroth order term of the perturbative development which is given by: 



t;^(o) 



-rrff 21, 



R,Z',MPo) = i^^ySip-Po)V - ipo) 



(84) 



and, V is the scattering matrix for a planar surface located at the height z — —H 



r^H21 



V {u) = exp{2iai{u)H)D2i{Po)- -D2i(Po 



^2i[Po)-\ ai(po)±a2(po) 



(85) 
(86) 
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The term exp{2 i ai{u) H) comes from the phase shift induced by the translation of the planar surface 
from the height z = to z — —H (see Appendix B2). The diagrammatic representation of Eq. ( |8^ ) is 
shown in Fig. |^, it is in fact a geometric series which can be summed 



U^'\p\p,) = {27:)H{p-p,)u'^'\p,) 



'^ -^-pW 



r=-TH21, 



U^°\p,)^ l + V--iPo)-V-'{Po) 



(87) 



From the previous results, Eq. (80) can be written in the following form : 

{0) —(0) —1 .— H 



U = U 



u^"^-v'-ar:,^^,^-u. 



where 












(89) 
(90) 



In order to obtain the perturbative development of Rd, we introduce the expansion 



R 



■H 



in (82) and 



which gives the following terms 



-'*-sei,e2 ~^ sei,e2 "'" -'•■sei,e2 



R 



H{3) 

'SCl,«2 ' 



(91) 



fl^ =F" + T^^f,;"^ !/<"'. t", 



Rd 



¥.u^^.R^:^^.u^^.¥^ 



R 



<2) 



t'-u^"^ 









fH{2) 



p".^. , R^(2) -(0) ^1 ^H(l) 
-"■s £i,£2 "■" -""s ei,e2 -''■sei,e2 



+ii— .^^"'.y^sf^;;^t7^"^F^sf^S;i 



■-^S£i,e2 ■ ^ 



(0) 



y fi 



if (2) 



.C7^''^T°, 



Using the development ( p3| ) for ii^^^^^, and the definitions (|70|-[7ll), (|7 
lations a development of the form (p3) for Rd with the following coefficients: 



(92) 
(93) 
(94) 



(95) 
we obtain after some calcu- 



-rrlO 



-— ff 21 



'zr ^rrlO / 



-rrff 21, 



^^(Po) - ( ^^"(Po) + V' '^(Po) ) • 1 + ^^"(Po) • ^" '^(Po) 



(96) 



this matrix is naturally diagonal, and its coefficients are identical to those of the reflection coefficients 
for a planar slabQ. The other coefficients are 



x['\p\p,)^t'\p)-u^'\p)-x'^11{p\p,)-u^"\p,)-t'\p,) , 
^f (pIpiIpo) -^'"(p) • t/^"^(p) • K;1(pIpiIpo) 



(97) 



-ai(pi)^fe!!e2(p|Pi) • U^"\pi) ■ V'\p,) ■ xfl]lipM)] ■ U^°\po) ■ t'\po) , 



^f (p|Pi|P2|Po) = t'\p) ■ U^"\p) ■ [xfl%{p\p,\pM) 

~ «l(P2)xf,;'^^(p|Pl|p2) • C/^"^(P2) •^'"(P2) •^f.!'l(P2|Po) 
-ai(Pl)^sei,e2(PlPl)-'^ (Pl)-y (Pl)-^sei,e2(PllP2lPo) 
+ ai(Pi)ai(p2)^sei,e,(p|Pl)-C^ (Pl)-'V^ (Pl)-^.ei.e2(Pl|P2) 
•f7'°\p2) •T^"(P2) • Xf,;i(p,|po)l . f/^^Po) .T-(p„) . 



(98) 



(99) 



In these expressions, X^^^ ^^(pIpq) = exp(i(ai(p) + ai(pQ)) i/) X^^^ g^(p|Po), and the subscripts ei, €2 



in X. 



means that we replace eo by ei, and ei by €2 in (pll 
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C A slab with a rough surface on the upper side 

We consider a slab delimited on the upper side by a two-dimensional rough surface, and on the bottom 
side by a planar surface, see Fig. m. To derive the reduced Rayleigh equation for this configuration, we 
have to combine the two following equations : 



<i^P -irri+'" 



i+,o- 



i2ny 



M," 'iu\p)-R^{p\p,)-E\p,) + M^"" {u\p,) ■ E\p,) = 



2(eo.i)^ar(«)^i+^^^ 



d'P inf^-o+ 



(2^) 



2^, 



iu\p) ■ Ru{p\Po) ■ E\po) + M„ ■ -{u\p,) ■ E\p, 



(ei - eo) 

2 {eg ei)^ ai{u) 
(ei - eo) 



(100) 



E^-{u). 



(101) 



with 



E^+{u)^v"^\u)-E^-{u), 



(102) 



-rH21 



where V is given by (|8^) , and i?„ is the global scattering matrix for the upper rough surface and the 
bottom planar surface. 

The reduced Rayleigh equation for the scattering matrix /?„ is then 



d^p 

(2^ 



M 



i+,o+ 

h 



-H21 



1-.0+, 



M. 



iu\p) + V {u)-M^' ^iu\p) ■ RuiplPo) = 



'{u\po) + V {u)-M^ ' (w|po) 



(103) 



With the expansion of I{a\p) in power series, we obtain the perturbative development: 

'1 



^u (Po) 



ai(Po)-"o(Po) ° 

i+,o- 



M 



(PoIPo) 



ai(Po) +"o(Po) 



V"'\P.) 



ai(Po) +ao(Po) 
-ai(po) +ao(Po) 



=-7lO, , -ttH21, ,\ r- T7IO/ X TTff21, , 

y {pg) + v (Po) • 1 + ^ (Po)-^ (Po) 



(104) 



^« (M|Po) = 2iQ (m|po) 
^u ("IPiIPo) = "i(«)Q (m|Po) + "o(Po)Q (m|Po) -2P (m|Pi)-Q (PiIPo) 
^« (w|Pi|P2|Po) = -3 [("?(«) + ao(Po))Q (M|Po) + 2ai(M)ao(Po)Q (w|Po) 



(105) 
(106) 



iP {u\p^)-X (pi|p2|Po) + * ai(M)P (M|p2)-ao(P2)-P (w|P2) -Q (P2IP0) 

(107) 



with 



Q °(w|Po) = 



2ao(Po) 



M mm) T7H21, , M (mm 

V (m) • — -— -— 

ai{u) + ao(wj 

i+,o- 



q;i(m) — q:o(m) 



aM (m|po)-X,,„_,^(po) 



M 



(w|Po) 



P^{u\p,)^ 



+6V (M).(aM (m|po).x1,;,,^(po) + M {u\p^^ 

m'+^°+(m|m) _h2i, , m'"'°+(m|m) 
V (m) 



ai(M) - ao(M) 

i+,o+ 



ai(M) + q:o(w) 



• M"""(M|pi)±y^''(M)-M' '"^(m|Pi) 



(108) 



(109) 
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where, a = zL, b — zL are the sign indices. After some computations we obtain: 

g++(w|po)-(ei-eo)[:D+(M)]-i- 

^ ei\\u\\\\po\\F+{u)F+{p„) -4Koa,{u)F-{u)F+{p,)iuxp,)A 

~eaai{u)ai{p„)Fy{u)Fy{pf^)u-pQ • [^"^ (p )]"S 

\-elKoa,{po)F+{u)Fy{po){uxpo), Kl F+{u) F+{p^)u ■ p^ j 



(110) 



Q (w|Po) = (ei -eo)[-Dio(M)]^ 



/ 



e,\\u\\\\p^\\Fy{u)F+{p^) 



-Co ai{u) ai(po) ^y (") Fy (Po) " ' Po 

1 

\-e§ Koai{p^)F^{u)F-{p^) {u x p„). 



i§ Koai{u)F+{u)F+{p„){u x p„),^ 
KiF^iu)F+ip„)u-p„ J 



[Dto{Po)r\ 



(111) 



Q^ (m|po) 



(^^-^°)[^ro(«r^ 



ao(Po) 

/ eoai{p„)\\u\\\\p„WF+{u)F-{p„) -e§ Koaiiu) aiip^) Fy {u)Fjj{pq) {ii x Pq)A 

-ei ai(M) a2(p^) Fy{u)F+{p„) u ■ p„ 

\-e„^ ei Koalipo) F+{u)F+{p„) {u x p^). 



Klai{p^)F+{u)F^{pQ)u-pQ j 



[DwiPoT' 
(112) 



Q (m|po) 



(ei "Co) r-pr+ , s.-l 



ao(Po) 

)IHIII 
-ei ai(M) a2(po) F+(M)i^+(po) « • Po 

V-eo"' ei ifoag(po)^H(M)^v(Po) (« x Po) 



/ eoai(Po)l|w|ll|Po!l^yN^y(Po) -^o ^o "!(") ai(Po) ^y (w)^H (Po) (" x Po)^^ 

-fCo "1 (Po) Fh iu)FH iPa)u-Pa J 



where 



and 



'^>^.4o))^(^-^^^^(-) 



T+V'\p,)-V"'\p,) 



[dUpoT' 

(113) 
(114) 



P^iu\p,)^ie,-eo)[Dl{u) 



P iu\p,) = (ei - eo) [£>io(w) 



/ \\u\\\\p,\\F+{u) 

+ai{u) aa {p^ ) Fy (u) u ■ p^ 

\eo ^ Xo ao(Pi) ^h(m) (« x Pi)^ 

/ ||m||||p||F^(m) 

+ai{u) ao(Pi) -Fy (m) m • Pi 

\eo ^ Xoao(Pi)^H(M)(M xpi). 



Elf ii'oal(M)F^^(M)(■u x Pi)^^ 

KlF+{u)u-p^ j 

(115) 

eJii'oa:i(M)i^^(M)(M X Pi)^^ 



if2J^i7N«-Pi / 



(116) 
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The first-order term was recently derived by Fuks et at3. Tliey have noticed that for this order, the 
matrix differs from the one obtained for a surface s eparating two semi-infinite media by only the factors 



ml) differ from Eq. 



F^. Lik ewise, fo r higher order, we see that Eqs. ( |110| 

for Eqs. ( |TT^ , [lT^ ) with respect to Eq. @, and Eqs. ( pl5[pl6| ) with respect to Eq. (|6^) 

thickness H becomes infinite, and the absorption Im(ei) ^ 0, or if ei = 62 

F — 1, and in that case we recover the matrix ( |66[]6q ) for a rough surface between two semi- infinite 

media. 



6|) by only _F , similarly 
So when the 

we have V =0, thus 



V. THE MUELLER MATRIX CROSS SECTION AND THE SURFACE STATISTIC 

When we consider an observation point in the far field iimit, the saddle-point method gives an asymp- 
totic form for the scattered field E" = E^^ obtained from Eq. ( pl| ) : 



^-(.,z)^ -p(:S,"-") 7(PiPo)-^'(Po), 



(117) 



with 



/(pIpo) = \ ■ -R(pIpo) : 

ZTTl 
P = Kq -TT-r: ' 



(118) 
(119) 



where 6 is the angle between e^ and the scattering direction (see Fig. g). In order to describe the incident 
and the scattered waves, we introduce the modified Stokes parameters: 



I'iP) 



I |i?f(p)P 

\EhkP)? 
2ReiEt.ip)Ef,ip)) 

V 2Im{E^yip)Ef,{p)) 



( 



I\P,) 



\ 



\EUPor 

\EhiP,)? 
2ReiElipo)E},ip,)) 

\ 2Im{E\.{p^)EUp^)) j 



(120) 



The analog of the scattering matrix for these parameters is the Mueller matrix, definedcl by: 

1 



1'{P) 



rMip\po)-r{Po), 



(121) 



which can be expressed as a functionEI of /(p|Po)- To maintain a matrix formulation in the following 
calculations, we introduce a new product between two-dimensional matrices with the definition : 



\JHV JHHJ \9hv 9hh 

/ fvvg'vv fvH9vH Re{fvv9vH) -Im{fvv9vH) 

^ fHv9Hv fHH9*HH Re{fHv9HH) -Irn{fHv9vH) 

I 2Re{fvv9*Hv) '^ReifyHghn) Re{fvvgvv + fHv9vH) -■^"^(/vvgW - Ivh9hv) 

\ 2Im{fvvg*Hv) '2ImifvH9HH) I'm-i.fvv9vv + fHvgvn) Re.{fvvgHH ~ fvHgnv) 

This product allows to express the matrix M as : 



(122) 



M(p|po) = /(p|po)0/(p|Po), 

K^ cos^ —, , , — , , , 

= (^^)2 Rip\Po) Q Rip\Po) ■ 



Following Ishimaru et al.cB, we define the Mueller matrix cross section per unit area a — (cry) 






(123) 
(124) 

(125) 





/ 1 1 1 


1 


110 1 


4 







\2 2 2 
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and, the bistatic Mueller matrixQ 7 = (7ij): 

7=- -M. (126) 

A cos 6*0 

These matrix are the generalization of the classical coefficients. In fact, if we assume, for example, that 
the incident wave is vertically polarized we have: 

J^I^Mp)P = j^liiiPlPoWviPr , (127) 

-^—^\Ef,{p)\' = 1 ^2l(p|Po)|S^(p)P • (128) 

A cos Wo ||r|p 

Thus 711 and 721, are respectively the classical bistatic coefficients jvv and jhv- We can also define the 
cross section and the bistatic coefficients for an incident circular polarization. As an example, taking the 
incident wave right circularly polarized, we have 

r(Po) = ^(l 10-2)*, (129) 

now, if we put a right-hand side polarizer at the receiver: 

_ _ , (130) 

\2 2 2/ 

the right to right bistatic coefficient jrr is: 

7rr = T (711 +712 + 2 714 + 721 +722+2 724 + 741 +742+2 744) , (131) 

where "fij are coefficients of the matrix 7 = (7^ ) . 

In a similar way we obtain the right to left bistatic coefficient: 

Jlr = T (711 +712+2 714 + 721 +722+2 724 - 741 "742-2 744) . (132) 

Up to now, we have made no hypothesis on the nature of the rough surface. Let us introduce the 
statistical caracteristics of the function h{x). We suppose that it is a stationary, isotropic Gaussian 
random process defined by the moments 

< h{x) > = 0, (133) 

<h{x)h{x')>^W{x~x'), (134) 

where the angle brackets denote an average over the ensemble of realizations of the function h{x). In 
this work we will use a Gaussian form for the surface-height correlation function W{x): 

W{x) = a^ exp{~x^/P) , (135) 

where a is the rms height of the surface, and I is the transverse correlation length. In momentum space 
we have: 

< h{p) > = 0, (136) 

< Hp) hip') > = (2^)2 Sip + p') Wip) , (137) 

with 

Wip) = I (fx Wix) expi-ip ■ x) (138) 

= Tr(j^l^ expi-p^l^/4). (139) 
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We are now able to define the bistatic coherent matrix 
1 



^coh _ 



A COS ^0 
and the incoherent bistatic matrix 



< /(pIpo) > < /(pIpo) >= . "^, _, < R{p\p,) > < R{p\Po) > , (140) 



yl(27r)2 cos6io 



r"'"^{p\Po) ^ 



1 



A cos 00 



[< f{p\Po) /(pIPo) > - < fiP\Po) > < fiP\Po) >] > 



^2 j,Qg2 ^ 

° [< fi(p|Po) -R(pIPo) > - < -R(pIPo) > < «(p|Po) >] • (141) 



A(27r)2 cos 6*0 
From Eq. (p3) , and the property of the Gaussian random process, we obtain 

rnh / I \ -/Vfi COS tf ^, . -=rC0h , , -=zCoh , , 

Y^'^iPlPo) = n Sip - Po) R (po) R (po) 

cos Uq 

R'"\po) ^ x'^"\po) + Ko cos 00 I^X^'\po\pM) W^(Pi - Po) 



(142) 
(143) 



:rC0h , 



where R (po) is a diagonal matrix describing the reflection coefficients of the coherent waves. For the 
incoherent part we have : 



7™-'^(p|Po) 



iiTo COS^ cos 00 r-=(l-l), I , t(2-2), |_ ^ , t(3-1), 



(2,r)2 



[/- ^^(p|Po) + r -(p|po) + r^ "(PlPo)] 



where 

y(l-l) 

t(2-2) 



(p|Po)^l^(p-Po)X^'^(p|Po)©X^'^(p|Po) 



(144) 



(145) 



(27r) 



(pIPo)^ / ^W^(p-Pi)mPi-Po)^''^(p|Pi|Po)0 [^''^(p|Pi|Po)+^*'\p|P+Po-PilPo) 

(146) 

(147) 



7^'-'^(p|Po) ^ W[p po) fx^'^(plpo) X^'^(pIpo) + x'''(p|po) X^'^(p|po)' 



with 



^^'\p|Po) 



(2^)2 



(3) 

W^(Pi-Po)^ (pIPoIPiIPo) 
+W^(p-Pi) fx'"(p|Pi|Po-p + Pi|Po)+^^'\p|Pi|p|Po) 



(148) 



VI. APPLICATIONS 



In the previous sections we have developed a method to compute the scattering matrices for a rough 
surface between two media, and for a thin film which includes one rough su rface. In this section, we 
will evaluate numerically the incoherent bistatic coefficients given by (144-147) for different values of the 
parameters which characterize the configurations. In all numerical simulations the media will be the 
vacuum (eo — 1). 

A rough surface separating to different media 

We consider that a polarized light of wavelength A = 457.9 nm is normally incident (0o = 0°, 0o — 0°) on 
a two-dimensional rough silver surface (see Fig. 0) characterized by the roughness parameters a — h nm, 
I = 100 nm, e = — 7.5 + ^0.24. As a matter of comparison we have chosen the same parameters used for the 
scattering by a one-dimensional rough surfaceEil. The perturbative development is given by Eqs. ( |6l| , |6^ . 
In Fig. ^, we present the results for an incident wave linearly polarized, the scattered field being observed 

in the incident plane (0 = 0°). The single scattering contribution associated with the term / is 
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(2 — 2) 

plotted as a dotted line, the double-scattering contribution I as a dashed line, the scattering term 

I as a dash-dotted line, and the sum of all these terms ^™'=°'' by the solid curve. 

We observed an enhancement of the backscattering which corresponds to the physical process in which 
the incident light excites a surface electromagnetic wave. In fact, the surface polariton propagates along 
the rough surface, then it is scattered into a volume wave due to the roughness, at the same time, a 
reverse partner exists with a path travelling in the opposite dire pt ii on , . These two paths can interfpe 
constructively near the backscattering direction to produce a peakllZlll§ll3. However, in one dimensioned, 
this peak can only be observed for a (TM) polarized incident wave because a surface polariton only exists 
for this polarization. In two dimensions, the surface wave also exists for a (TM) polarization, but in fact 
a depolarization occurs so that a (TE) incident wave can excite a (TM) surface wave, and this surface 
wave can be scattered into volume wave with both polarizations as can be seen in Fig. pi Now, when the 
incident wave is circularly polarized, we see in Fig. |9|, that the enhanced backscattering takes also place. 
We have not displayed the left to left, and left to right polarizations because the media are not optically 
active, as a consequence, the results are the same either the incident wave is right or left polarized. 

r(2-2) 



In the expression (144), the peak is produced by the term I . We see that the term P{u\pi) 



Q^{Pi\Po) i^ ^sea.eiiPlPilPa) Contains a factor of the form (see Eq. (|5^)): 



1 

eiQ:o(Pi) + eoai(Pi)' 

which is close to zero excepted when Pi is near the resonance mode p^ of the polariton, which is given by 
)+ ( 



[DtovviPi)r' = —r7:r^—r7zrT^ (149) 

3n Pi i; 
the roots D^Qy{pj.) = 0. When we observe the field scattered far away from the backscattering direction 



f2) (2) 

(P+Po + 0)> the terms X,^^^^^(p|Pi|Po) and X^^^ ,,^(p|p-HPo-Pi|Po) conteningD+ ^ are non zero when 
P\ ~Pr^ ^-iid P + Po ~ Pi ~Pr respectively. Since these domains are disjoint, the product of these two 
terms is approximatively zero. Conversely, when we are near the backscattering direction (p -f Pg ~ 0), 
the terms inside brackets are almost equal and produce the enhancement factor. This enhancement 

factor is not equal to 2 because the matrices Q (p|Po) and Q (p|Po) in -X'^g^^ ^^ do not contain the term 
[I?j*Qy(P]^)]~"'^, so they produce a significative contribution whatever the scattering angle is. In order to 
isolate more precisely the terms producing an enhanced backscattering, a better approach is to work with 
the formalism of Re f. [iq derived from quantum mechanical scattering theory, such an approach is used 
for instance in Ref. [L7|J23[ If the decomposition of each step of the mutiple scattering process is clearly 
put in evidence, however, it offers the disadvantage to produce a more heavier perturbative development 
as it can be seen when comparing (^) and ( p6[J6q ) with (15-19, A-1) of Ref. ^. 

A film with a rough surface on the upper side 

We consider a dielectric film (see Fig. 0) of mean thickness H — 500 nm, dielectric constant 
t\ — 2.6896 -I- 0.0075i, deposited on a planar perfectly conducting substrate (ti = — oo) and illumi- 
nated by a linearly polarized light of wavelength A = 632.8 nm normaly incident {(/jq = 0°, Oq = 0°). The 
two-dimensional upper rough surface is characterized by the parameters a = 15 nm and I = 100 nm. The 
scattering diagrams are shown in Fig. [l O| with the sarne cu rve labelling as before. The perturbative de- 
velopment being given by Eqs. (|l05| , 107| ) a nd Eq s. ( 110| , 116|) . Since we have chosen an infinite conducting 



plane (e2 — — cxd) the coefficients F^ (Eq. (114)) have the following form 



-.±f x_ l±exp(2iao(Po)-ff) 



^ ° (eiao(Po)+eoai(Po)) + (eiao(Po)-eoai(Po)) exp(2iao(Po)-ff) ' 

F±(„)^ lTexp(2iao(po)g) 

^^■^°^ (ao(Po) +ai(po)) + (ao(Po) " «i(Po)) exp(2iao(Po) ^) ' 



The parameters are the same as those used in Ref. 23 for a one-dimensional dielectric film where a (TE) 



polarized wave is incident. The thickness was chosen in way that the slab supports only two guided wave 
modes, p}p^ = 1.5466 iiTg, and p'^^ = 1.2423 -ftTo for the {TE) polarization. These modes are resonance 
modes, they verify [F^]~^{prp^) — 0. For the (TM) case, we have three modes given by the roots 

[F^r^iETM) = 0, which are : p^j^/ = 1-6126 i^o, Ptm = 1-3823 i^o and p^j^j = 1.0030/^0- As described 
in Ref. [lipi^p7| , these guided modes can produce a classical enhanced backscattering with satellite peaks 
symmetrically positioned. The satellite peaks angles are given by the equation : 

sin ^r ^ - sin 00 ± 1^ b" - Pn , (152) 



Ko 
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where p", p™ describe one of the guided mode, when n = m, we recover the classical enhanced backscat- 



tering. We can give an explanation of this formula as in the previous case. In the expression (106), the 
term producing the peaks comes from 

-2P(m|pi)-Q^(PiIPo) (153) 

where Q (Pi|Po) contains the factors F^{pi) having resonances for the slab guided mode. The product 

^f (p|Pi|Po)0^f (pIp+Po -PilPo) (154) 



in Eq. (146) has a significant contribution only when Pi and p + Po ^ Pi ^^^ near resonance modes. 
As there are several resonances, we can have Pi w p^ and p + Po ^ Pi ~ Pm with n ^ m where p„ 
and p„ are resonances vector. If p„ = zLpne^ and p„ = TPine^ (guided modes propagating along 
the incident plane but with opposite directions), we have (p + pp) • e^. w ±(Pn — Pm) which is another 



way of writing Eq. (152). For the {TE) polarization, since we have only two guided waves, the satellite 
peaks can only exist at the angles 0]^{TE) = ± 17.7°. Now, for the (TM) polarization we have three 
possibilities: e^^{TM) = ±13.3°, e^^{TM) = ±37.6°, and ef{TM) = ±22.3°. The satellite peaks are 

(2—2) 

produced by the term I , in the case of (TM) polarization we do not get any significant contribution 
to satellite peaks. However, for the (TE) to [TE) scattering shown in Fig. O, we find satellite peaks 
at the angle 9}^(TE) — ± 17.7° positioned along a dotted line. Now, by doubling the slab thickness, see 
Fig. O, the satellitajDeaks disappear for all the polarization, but we see a new phenomenon called the 
Selenyi fringes .tZfotj. For a slightly random rough surface, the slab produces fringes similar to those 
obtain with a Fabry-Perrot interferometer illuminated by an extended source. The roughness modulates 
amplitude fringes but their localization remains the same as for the interferometer. We also notice that 
the enhanced backscattering decreases with the slab thickness. We can conclude that as in the case of 
one-dimensional rough surface, the satellite peaks appear only when the wave guide supports few modes 
for the (TE) polarizations. These results differ from those obtained in Ref. ^ where no satellite peak 
appears in their two-dimensional slab, we have checked that with their parameters values we also find 
no peak, and we agree with the results given by the contributions of the first and second order terms. 
However, the third order gives a contribution larger than the first one, such a result casts some doubt on 
the validity of the (SPM) method in that case. 

However, for the choice of parameters presented here no satellite peak has been observed even when the 
thickness of the slab is chosen in such a way that only two guided modes exist for the (TM) polarization 
(a result not presented here). This is in agreement with the results of Ref. B^ for one dimensional surface 
where it is noticed that the excitation of [TM) modes are more difficult to excite than the (TE) modes. 
In order to enhance this effect they choose a higher permittivity for the media 1 : ei = 5.6644 + i0.005. 
In this case satellite peaks are observed for a slab which supports three guided modes. We have also 
done numerical calculations with these parameters, however we do not observe satellite peaks. So the 
transition from one dimensional to two dimensional rough surface lower the efficiency of the excitation 
of (TM) modes. Next, instead of doubling the slab thickness, we have changed the infinite conducting 
plane by a silver plane (e2 = —18.3 + 0.55i), we see in Fig. O, that the enhancement of backscattering is 
also decreased, and that there is no more satellite peak corresponding to (TE) to (TE) scattering. This 
fact has to be compared with the next configuration, where the rough surface is now between the media 
1 and the media 2, see Fig. ^. 

A film with a rough surface on the bottom side 

The permitivities are the same as in the previous configuration, excepted that the case £2 — —00 cannot 

be treated with the (SPM) because the second and third order diverge. The rms height a has now the 

value cr = 5 nm and I = 100 nm. We have not chosen cr = 15 nm because numerically we have noticed that 

—(1-1) —(2-2) 

the first order term I was not greater than the second order I , which means that we are near 

the limit of validity of (SPM). The perturbative development is given by (|97| - P9|) , and the guided modes 
are the roots of [Xyy j^{ptm)]^^ for (TM) polarization, and of [Xj^lj j^{pte)]~^ for (TE) polarization. 
We obtain two modes in the (TE) case, whose values are : p}p^ — 1.5534^0, and p'^^ — 1.2727 A'o, 
the corresponding satellite peaks angles are : 9}^{TE) — ±16.3°. For the [TM) case, we have three 
guided modes with p^j^.j = 1.7752 i^Tp, p^j^i = 1.4577 7^0 and p^j^j — 1.034 ii'o, they correspond to six 
possible sattelite peaks angles given by e^^{TM) = ±18.51°, e^^{TM) = ±47.8°, and 0f{TM) ^ ±25°. 
We see the apparition of satellite peaks only for (TM) to (TM) scattering process as shown in Fig. ^. 



19 

This result differs from the previous case because, on one hand, the rough surface being not a perfect 
conductor we still obtain satellite peaks, on the other hand, these satellite peaks now appear for the (TM) 
to (TM) polarization instead of (TE) to (TE). This is a surprising result because the {TM) polarization 
which has one more mode than the {TE) one, should decrease the amplitude of the satellite peaks for 
this polarization as it was the case with an upper rough boundary. Moreover, we see in Fig. |l^, that the 
three satellite peaks can be clearly separated. This can be explained from the fact that there are two 
phenomenons which occur in this case. The first is the same as in the previous case where the wave can 
excite guided modes through the roughness which produces the enhancement of backscattering and the 
satellite peaks. These effects come from the term 

a,{p{)X",^'l{p\p,)-U^'\p,)-v'\p,)-x"}'l{p,\p,) (155) 

in Eq. ( JQ^ ) where U {pi) have resonances for the different modes of the guided wave. But, there is 
also a second phenomenon which was described in our first example where the rough surface can excite 



a plasmon mode. This appears from Eq. (98) with the term 



^f.il(PIPilPo) (156) 



and subsequently in Eqs. (144,146). The localization of this mode p^ is given by : 



[DtivviPr)]'^ = 7-T^ 1-^ ■ (157) 

L 2ivvyi^ri\ €2ai{p,)+eia2{Pr) 

In o ur ca se this give Hp^H — 1. 7755 - RTp w hich is very close to the value p^jv/ — 1-77^2 Kq. So, in 
Eq. ( |146| ) the product of (155) by (156) can produce peaks where p^ ~ ±||p^||ej; « ±^^^,^63, and 



(Po +P) ■ ^x = i(PTfl/ ~Ptm) with n — 1,2, 3. We have effectively verified numerically that the product 
of this two term can enhanced considerably the different peaks in particular the first satellite peaks 6}^ 
(when n — 2). Now, by doubling the slab thickness, we see from Fig. |l^ that the satellite peaks have 
disappeared due to the too many guided modes which can be excited. 

VII. CONCLUSIONS 

We have obtained four generalized reduced Rayleigh equations which are exact integral equations, and 
where one of the four unknown fields coming on the rough surface has been eliminated. These equations 
offer a systematic method to compute the small perturbation development without lenghty calculations, 
moreover, the scattering matrices are only two dimensional. All the theoretical calculations have been 
done up to order three in the height elevation which allow us to obtain all the fourth-order cross-section 
terms. We have calculated the perturbative development for three different structures composed of a rough 
surface separating to semi-infinite media, and a dielectric film where one of the two boundaries is a rough 
surface. For the first structure, the perturbative expression has been already calculated at the third order, 
but our derivation offers the advantage to be formulated in compact manner making easier numerical 
computations. For the slabs configuration we present new results. It has to be noticed that for the case 
of a rough surface in the upper position, the generalized derivation of the reduced Rayleigh equations 
becomes mandatory. The numerical results show an enhancement of the backscattering for co- and cross- 
polarizations in all these cases. In the slab case, for some configurations and definite polarizations, we 
have detected satellite peaks which result from interference of different waveguide modes. This general 
formulation can be extended to the configuration including two rough surfaces, and some results will be 
presented in a next paper. 
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APPENDIX A: THE INTEGRATION BY PARTS 



We need to calculate the following integral: 



(Al) 



Since Vh{x,y) is zero for \x\ > L/2 or \y\ > L/2, we can fix the integration limits. We choose the 
boundary limits xi in x such that |a;;| > L/2, and {u — p)x xi = 2'Kmx, with tUx € Z. Similarly, we choose 
the boundary yi in y such that \yi\ > L/2 and {u — p)yyi — 2'Kmy, with my £ Z. Thus, the integral ( |A1[ ) 

is : 

/ dxdy exp{—i{u — p)-x) Vh{x) exp{~i{bai{u) ~ aai{p))h{x)) 

-x, J -Vi 

exp(— i(M — p) ■ X — i{bai{u) — aai{p))h{x))'^^~ ^' 



/Vl 
dy 

/xi 
dx 
-xi 
xi rVi 



~i{bai{u) — aai{p)) 
exp{—i{u — p) ■ X — i{bai{u) — aai{p))h{x)) 



-i{bai{u) — aai{p)) 



x=-xi 

y=+yi 



-i(M-p) 



-xiJ-yi iibai{u) - aaiip)) 
2™ (u-p) 



exp{—i{u — p) ■ X — i{bai{u) — aai{p))h{x)) 



= - d'x 



(bai{u) — aai{p)) 



eM-^ikl' kl^) ■ r.) , 



(A2) 
(A3) 



The term in the square bracket canceled due to the choice made for xi and yi. From the previous 
calculations, we can now replace \Jh{x) by: 



Vh{x) 



{u-p) 



{bai{u) — aai{p)) 



(A4) 



APPENDIX B: PERTURBATIVE DEVELOPMENT AND RECIPROCITY CONDITION 

As was noticed by Voronovich,a the scattering operator R has a very simple law of transformation 
when we shift the boundary in the horizontal direction by a vector d: 



Rx^h(a^-d){p\Pa) = cxp[-i(p - Po) • d] R^^h{j:){p\Po) -. 
or, when we translate the surface by a vertical shift H e^' 

Rh+Hip\Po) = exp[-i(ao(p) + ao(Po)) H] Rh ■{p\Pq) 



(Bl) 



(B2) 



Now, using (Bl), we can deduce some properties on the perturbative development of the scattering 
operator. The generalization of the Taylor expansion for a function depending on a real variable to an 
expansion depending on a function (which is in fact a functionnal) can be expressed in the following form: 



R{p\p,) = R^'\p\po) + R^'\p\Po) + R^'\p\Po) + R^'\p\Po) 



where 



R^'\p\Po)=I^R^'\p\pM)h{p^), 

R {p\Pi\P2\Pa) HPi) HP2) . 

R iP\Pi\P2\P3\Po) HPi) hiP2) HPs) 



r'''\p\po) = 



R {p\Po 



(27r)2 

dVi d^P2 d(2) 
(27r)2 (27r)2 

d^Pi d2p2 d2p3 (3) 



(27r)2 (27r)2 {2tt) 



(B3) 

(B4) 
(B5) 
(B6) 
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Applying this perturbative development on each side of (Bl), and taking their functional derivative (see 
Rcf. 0) defined by: 



Sin) 



5h{q{)..5h{q^) 



We obtain for all ?i > in the limit /i = 

TJ(") 



9„ - Po) ■ <^) -R " (pIqiI 



R (p|gil---|g„|Po) = exp(-i(p-qi 

We find that 

R (p|Qi|---|g„|Po)oc(5(p-qi 
so we can define X matrices by the relations : 

R {p\Pq) ^ otQ{pfi) I -—^ X {p\Pi\Po)h{p~p^)h{p^-pfi), 

d^Pi d2p2 ^(3) 



-Po) 



knlPo) 



(B7) 

(B8) 
(B9) 

(BIO) 
(Bll) 

(B12) 



r'-'\p\p,) - ao(Po) // ^^ ^^'^(pIPiIp^IPo) Mp " Pi) MPi - P2) MP2 " Po) , (B13) 



where ao(Po) i^ introduced for a matter of convenience. 

Let us now make some remarks about the reciprocity condition. If we define the anti-transpose opera- 
tion by: 



a b 
c d 



aT 



a —c 

-b d 



the reciprocity condition for an incident and a scattered waves in the medium readgj: 



■^aT 



R (pIpo) _ fl(-Po|-p) 

ao(Po) 



ao(p) 



(B14) 



(B15) 



Making use of the previous functional derivative, we would like to prove that each order of the perturbative 
development must satisfies this condition. It is easy to show that 



"-^(1) 



>|Po) 



^^'"(-Pol-P), 



(B16) 



thus -X" is reciprocal, but the same conclusion cannot be extended to X when n > 2. For example. 



in the case n = 2, using ( B15| ), we can only deduce that: 

(27r)2 



^ (pIPiIPo) Hp - Pi) HPi - Po) 



d^Pi ^(2 



{2n) 



From this we cannot deduce a result similar to (B16) for X . This fact is well illustrated with the 
following identity(which can be demonstrate with a transformation of the integration variables)) : 



2^ (-P0I-P1I -p)^(p-Pi)/^(Pi-Po)' 

(B17) 



r(2) 



(2^)2 



{p+Po-'^Pi)Hp-Pi)KPi -Po) =0. 



(B18) 



We see that Pi — > p + Pq — 2pi is not the null function although the integral is null. From this we 



r(«) 



deduce that X for n > 1 are not unique. Moreover in using (B18) we can transform the X in a 



r(«) . 



reciprocal form. This procedure is illustrated in the one-dimensional case in Ref. ^, and the results for 
the second-order in the electromagnetic case are given in Ref. || and Ref. |2^. 
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FIG. 1: A rough surface with an incident wave coming from both sides of medium and 1. 
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FIG. 2: Decomposition of the wave vector k^ 
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FIG. 3: A two-dimensional rough surface separating two dielectric media and 1. 
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FIG. 4: A slab formed with a bottom two-dimensional rough surface and an upper planar surface. 
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FIG. 5: Definitions of scattering matrices for a planar surface. 
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FIG. 6: Diagrammatic representation of the operator U 
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FIG. 7: A slab formed with an upper two-dimensional rough surface and a bottom planar surface. 
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FIG. 9: The same configuration as Fig. 0, but with a right incident cicularly polarized wave, and a right 
to right (or left to left), right to left (or left to right) observed polarizations 
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FIG. 10: The bistatic coefficients for an liorizontal {TE) and a vertical (TM) polarized incident light of 
wavelength A — 632.8 nrri, on a slab with an upper two-dimensional randomly rough surface, characterized 
by the parameters a = 15 nm, I = 100 nm, ei = 2.6896 -I- i0.0075, thickness H = 500nm, deposited on 
an infinite conducting plane (e2 = — oo). The scattered field is observed in the incident plane. For each 



the total incoherent scattering 7 

(2 — 2) — 

(dotted line), the second order I (dashed line), and the third order / 



■mcon ^gQjjj curve), the first order given by / 

t(3-i) 



(dash-dotted line). 
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FIG. 11: Details of second order (TE) to (TE) contribution to the scattering shown in Fig. |T^. We see 
two sateUite peaks at the angle 9^{TE) = ± 17.7°, the dotted-lines mark the peaks angle position. 
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FIG. 12: Effect of the slab thickness, H — lOOOriTO, on the configuration shown in Fig. HC 
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FIG. 13: The same parameters as in Fig. |o[ but with a silver plane characterized by 62 = —18.3 + 0.55i. 
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FIG. 14: The bistatic coefficients for an liorizontal (T-E) and vertical (TM) polarized light of wavelength 
A = 632.8 n?Ti, incident on a film of permittivity ei — 2.6896 + i0.0075, deposited on a two-dimensional 
randomly rough surface, characterized by the parameters, a = 5nm, I = 100 nm, €2 = —18.3 + 0.55i, 
thickness H — 500nm. The scattered field is observed in the incident plane. For each figure are plotted 

: the total incoherent scattering ^*"'=°'' (solid curve), the first order given by I (dotted line), the 

(2 — 2) (3—1) 

second order I (dashed line), and the third order / (dash-dotted line). 
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FIG. 15: Details of the second order {TM) to (TM) contribution to the scattering shown in Fig. p^ , 
dotted Unes mark the peaks angle position. 
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FIG. 16: Effect of the slab thickness, H — lOOOriTO, on the configuration shown in Fig. UA 
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